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INTRODUCTION 



1.1 A CONTROL ENGINEERING APPROACH TO FUZZY CONTROL 

This book gives a comprehensive treatment of model-based fuzzy control 
systems. The central subject of this book is a systematic framework for the 
stability and design of nonlinear fuzzy control systems. Building on the 
so-called Takagi-Sugeno fuzzy model, a number of most important issues in 
fuzzy control systems are addressed. These include stability analysis, system- 
atic design procedures, incorporation of performance specifications, robust- 
ness, optimality, numerical implementations, and last but not the least, 
applications. 

The guiding philosophy of this book is to arrive at a middle ground 
between conventional fuzzy control practice and established rigor and sys- 
tematic synthesis of systems and control theory. The authors view this 
balanced approach as an attempt to blend the best of both worlds. On one 
hand, fuzzy logic provides a simple and straightforward way to decompose 
the task of modeling and control design into a group of local tasks, which 
tend to be easier to handle. In the end, fuzzy logic also provides the 
mechanism to blend these local tasks together to deliver the overall model 
and control design. On the other hand, advances in modern control have 
made available a large number of powerful design tools. This is especially 
true in the case of linear control designs. These tools for linear systems range 
from elegant state space optimal control to the more recent robust control 
paradigms. By employing the Takagi-Sugeno fuzzy model, which utilizes local 
linear system description for each rule, we devise a control methodology to 
fully take advantage of the advances of modern control theory. 
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INTRODUCTION 



We have witnessed rapidly growing interest in fuzzy control in recent 
years. This is largely sparked by the numerous successful applications fuzzy 
control has enjoyed. Despite the visible success, it has been made aware that 
many basic issues remain to be addressed. Among them, stability analysis, 
systematic design, and performance analysis, to name a few, are crucial to the 
validity and applicability of any control design methodology. This book is 
intended to address these issues in the framework of the Takagi-Sugeno 
fuzzy model and a controller structure devised in accordance with the fuzzy 
model. 



1.2 OUTLINE OF THIS BOOK 

This book is intended to be used either as a textbook or as a reference for 
control researchers and engineers. For the first objective, the book can be 
used as a graduate textbook or upper level undergraduate textbook. It is 
particularly rewarding that using the approaches presented in this book, a 
student just entering the field of control can solve a large class of problems 
that would normally require rather advanced training at the graduate level. 

This book is organized into 15 chapters. Figure 1.1 shows the relation 
among chapters in this book. For example, Chapters 1-3 provide the basis 
for Chapters 4-5. Chapters 1-3, 9, and 10 are necessary prerequisites to 




Research level 
(research topics) 



Fig. 1.1 Relation among chapters. 
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understand Chapter 11. Beyond Chapter 3, all chapters, with the exception of 
Chapters 7, 11, and 13, are designed to be basically independent of each 
other, to give the reader flexibility in progressing through the materials of 
this book. Chapters 1-3 contain the fundamental materials for later chapters. 
The level of mathematical sophistication and prior knowledge in control have 
been kept in an elementary context. This part is suitable as a starting point in 
a graduate-level course. Chapters 4-15 cover advanced analysis and design 
topics which may require a higher level of mathematical sophistication and 
advanced knowledge of control engineering. This part provides a wide range 
of advanced topics for a graduate-level course and more importantly some 
timely and powerful analysis and design techniques for researchers and 
engineers in systems and controls. 

Each chapter from 1 to 15 ends with a section of references which contain 
the most relevant literature for the specific topic of each chapter. To probe 
further into each topic, the readers are encouraged to consult with the listed 
references. 

In this book, S > 0 means that 5 is a positive definite matrix, S > T 
means that S — T > 0 and W = 0 means that W is a zero matrix, that is, its 
elements are all zero. 

To lighten the notation, this book employs several particular notions which 
are listed as follow: 



i < j s.t. hj n hj =£ </>, 
i < j s.t. h t n hj + 4>. 

For instance, the condition (2.31) in Chapter 2 has the notation, 

i <j < r s.t. h t n hj ¥= 4>. 

This means that the condition should be hold for all i < j excepting li, n hj 
= 4> [he., /i ; (z(t)) X hj(z(t)) = 0 for all z (?))], where /r,(z(t)) denotes the 
weight of the ith rule calculated from membership functions in the premise 
parts and r denotes the number of if-then rules. Note that h L n hj = <$> if 
and only if the ith rule and j th rule have no overlap. 
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TAKAGI-SUGENO FUZZY 
MODEL AND PARALLEL 
DISTRIBUTED COMPENSATION 



Recent years have witnessed rapidly growing popularity of fuzzy control 
systems in engineering applications. The numerous successful applications 
of fuzzy control have sparked a flurry of activities in the analysis and design 
of fuzzy control systems. In this book, we introduce a wide range of analysis 
and design tools for fuzzy control systems to assist control researchers and 
engineers to solve engineering problems. The toolkit developed in this book 
is based on the framework of the Takagi-Sugeno fuzzy model and the 
so-called parallel distributed compensation, a controller structure devised in 
accordance with the fuzzy model. This chapter introduces the basic concepts, 
analysis, and design procedures of this approach. 

This chapter starts with the introduction of the Takagi-Sugeno fuzzy 
model (T-S fuzzy model) followed by construction procedures of such models. 
Then a model-based fuzzy controller design utilizing the concept of “parallel 
distributed compensation” is described. The main idea of the controller 
design is to derive each control rule so as to compensate each rule of a fuzzy 
system. The design procedure is conceptually simple and natural. Moreover, 
it is shown in this chapter that the stability analysis and control design 
problems can be reduced to linear matrix inequality (LMI) problems. The 
design methodology is illustrated by application to the problem of balancing 
and swing-up of an inverted pendulum on a cart. 

The focus of this chapter is on the basic concept of techniques of stability 
analysis via LMIs [14, 15, 24]. The more advanced material on analysis and 
design involving LMIs will be given in Chapter 3. 
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6 TAKAGI-SUGENO FUZZY MODEL AND PARALLEL DISTRIBUTED COMPENSATION 

2.1 TAKAGI-SUGENO FUZZY MODEL 

The design procedure describing in this book begins with representing a 
given nonlinear plant by the so-called Takagi-Sugeno fuzzy model. The fuzzy 
model proposed by Takagi and Sugeno [7] is described by fuzzy IF-THEN 
rules which represent local linear input-output relations of a nonlinear 
system. The main feature of a Takagi-Sugeno fuzzy model is to express the 
local dynamics of each fuzzy implication (rule) by a linear system model. 
The overall fuzzy model of the system is achieved by fuzzy “blending” of the 
linear system models. In this book, the readers will find that many nonlinear 
dynamic systems can be represented by Takagi-Sugeno fuzzy models. In fact, 
it is proved that Takagi-Sugeno fuzzy models are universal approximators. 
The details will be discussed in Chapter 14. 

The ith rules of the T-S fuzzy models are of the following forms, where 
CFS and DFS denote the continuous fuzzy system and the discrete fuzzy 
system, respectively. 

Continuous Fuzzy System: CFS 

Model Rule i: 



IF z,(t) is A/,, and ••• and z p (t) is M ip , 

(x(t) =A/x(t) +B:u(t), 
THEN { , ' ' ) ' 

\y(t) = C t x(t), 

Discrete Fuzzy System: DFS 

Model Rule i: 



i = 1, 2, . . . , r. 



( 2 . 1 ) 



IF z,(t) is M n and ••• and z p (t) is M ip , 



THEN 



x(t + 1) = A i x(t) + BjU(t), 
y(t) = CjX(t), 



i = 1 , 2 ,..., 



( 2 . 2 ) 



Here, M,- ■ is the fuzzy set and r is the number of model rules; x(t) e R n is 
the state vector, u(t ) e R m is the input vector, y(t ) e R q is the output 
vector, A i <ER nXn , B i ^R nXm , and C i ^R qXn \ z/t), . . . , z p (t) are known 
premise variables that may be functions of the state variables, external 
disturbances, and/or time. We will use z(t ) to denote the vector containing 
all the individual elements z/t), . . . , z p (t). It is assumed in this book that the 
premise variables are not functions of the input variables u(t). This assump- 
tion is needed to avoid a complicated defuzzification process of fuzzy 
controllers [12]. Note that stability conditions derived in this book can be 
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applied even to the case that the premise variables are functions of the input 
variables u(t). Each linear consequent equation represented by Ajx(t) + 
Bjiiit) is called a “subsystem.” 

Given a pair of (jc(t), u(t)), the final outputs of the fuzzy systems are 
inferred as follows: 



CFS 

X>.-(z(0)M»*(0 + B M 0} 

x(t) = — r 

I>;(z(0) 

i = 1 

= X>,-(z(0)M.-*(0 + B M { )}’ ( 2 - 3 ) 

i = 1 
r 

'Lw i {z(t))C i x{t) 

y(t ) = —r 

E w,-(z(0) 

i = 1 

= f>i(z(0) c »*(0- ( 2 - 4 ) 

/ = 1 

DFS 

r 

Hwi(z(t)){A iX (t) + Bjii(t)) 

X(t + 1) = — ■ 7 

Z>;(z(0) 

i = 1 

= I>;(z(0)K*(0 + B ««(0}> (2-5) 

;=1 

r 

Ew.-(z(0)C/*(0 

J(0 = iZL - F 

E^(z(0) 

/ = 1 

= EMz(0) c i*(0» 

( = t 



( 2 . 6 ) 
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where 



z(0 = [^l(0 z 2(0 - Zp(t)]> 

W/(z(0) = n^,( z ;(0)> 

7 = 1 

A i(z(0) = — 

£w.-(z(0) 



(2.7) 



for all t. The term M,--(z-(0) is the grade of membership of zff . ) in Af,--. 
Since 



we have 



| E w,-(z(0) > °. 

|w,(z(t)) >0, i = 1,2,. . . , r. 



( 2 . 8 ) 



for all t. 



| EMz(O) = i» 

|/i,(z(f)) > 0, 7 = 1,2, ...,r, 



(2.9) 



Example 1 Assume in the DFS that 
p = n, 

z x (t) = x(t), z 2 (t) = x(t — 1),. . ., z n (t) = x(t — n + 1). 
Then, the model rules can be represented as follows. 

Model Rule i: 



IF x(t) is M n and ••• and x(t — 77 + 1) is M in , 

l X(t + 1) = A,JC(t) + 

THEN ) , [ 7 = 1, 2, 

\j(0 = c,.*(0, 

where x(t) = [x(t) x(t — 1) ••• x(t - n + l)] r . 



Remark 1 The Takagi-Sugeno fuzzy model is sometimes referred as the 
Takagi-Sugeno-Kang fuzzy model (TSK fuzzy model) in the literature. In this 
book, the authors do not refer to (2.1) and (2.2) as the TSK fuzzy model. The 
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reason is that this type of fuzzy model was originally proposed by Takagi and 
Sugeno in [7]. Following that, Kang and Sugeno [8, 9] did excellent work on 
identification of the fuzzy model. From this historical background, we feel 
that (2.1) and (2.2) should be addressed as the Takagi-Sugeno fuzzy model. 
On the other hand, the excellent work on identification by Kang and Sugeno 
is best referred to as the Kang-Sugeno fuzzy modeling method. In this book 
the authors choose to distinguish between the Takagi-Sugeno fuzzy model 
and the Kang-Sugeno fuzzy modeling method. 



2.2 CONSTRUCTION OF FUZZY MODEL 

Figure 2.1 illustrates the model-based fuzzy control design approach dis- 
cussed in this book. To design a fuzzy controller, we need a Takagi-Sugeno 
fuzzy model for a nonlinear system. Therefore the construction of a fuzzy 
model represents an important and basic procedure in this approach. In this 
section we discuss the issue of how to construct such a fuzzy model. 

In general there are two approaches for constructing fuzzy models: 

1. Identification (fuzzy modeling) using input-output data and 

2. Derivation from given nonlinear system equations. 

There has been an extensive literature on fuzzy modeling using input-out- 
put data following Takagi’s, Sugeno’s, and Kang’s excellent work [8, 9]. The 
procedure mainly consists of two parts: structure identification and parame- 
ter identification. The identification approach to fuzzy modeling is suitable 




Fig. 2.1 Model-based fuzzy control design. 
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for plants that are unable or too difficult to be represented by analytical 
and/or physical models. On the other hand, nonlinear dynamic models for 
mechanical systems can be readily obtained by, for example, the Lagrange 
method and the Newton-Euler method. In such cases, the second approach, 
which derives a fuzzy model from given nonlinear dynamical models, is more 
appropriate. This section focuses on this second approach. This approach 
utilizes the idea of “sector nonlinearity,” “local approximation,” or a combi- 
nation of them to construct fuzzy models. 

2.2.1 Sector Nonlinearity 

The idea of using sector nonlinearity in fuzzy model construction first 
appeared in [10]. Sector nonlinearity is based on the following idea. Consider 
a simple nonlinear system x(t) = where /( 0) = 0. The aim is to find 

the global sector such that x(t) = f(x(t )) e [a : a 2 ]x(t). Figure 2.2 illustrates 
the sector nonlinearity approach. This approach guarantees an exact fuzzy 
model construction. However, it is sometimes difficult to find global sectors 
for general nonlinear systems. In this case, we can consider local sector 
nonlinearity. This is reasonable as variables of physical systems are always 
bounded. Figure 2.3 shows the local sector nonlinearity, where two lines 
become the local sectors under —d < x(t) < d. The fuzzy model exactly 
represents the nonlinear system in the “local” region, that is, —d < x(t) < d. 
The following two examples illustrate the concrete steps to construct fuzzy 
models. 



a t x(t) 




Fig. 2.2 Global sector nonlinearity. 
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Example 2 Consider the following nonlinear system: 

/ *i(0 ' 

\x 2 (t) / 



-Xi(t) +x 1 (t)x 3 2 (t ) 
-x 2 (t) + ( 3 + x 2 (t))x 3 (t) 



( 2 . 10 ) 



For simplicity, we assume that Xj(0 e [— 1,1] and x 2 (t) e [— 1 , 1] . Of 
course, we can assume any range for x^t) and x 2 (t) to construct a fuzzy 
model. 

Equation (2.10) can be written as 



x(t) 



-1 

(3 +x 2 (t))xf(t) 



Xi (t)xi(t) 
-1 



x(t), 



where jc(r) = [x^t) x 2 (t)] T and x x {t)x\{t) and (3 + x 2 (t))x 2 (t) are nonlinear 
terms. For the nonlinear terms, define Zj(r) = xfa)x 2 (t) and z 2 (t ) = (3 + 
x 2 (t))xj(t). Then, we have 



x(t) 



-1 

z 2 (t) 



Zl(t) 

-1 



x(t). 
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Next, calculate the minimum and maximum values of z A (t) and z 2 (t ) under 
Xj(t) e [ — 1,1] and x 2 (t ) e [-1, 1]. They are obtained as follows: 

max Zj(t) = 1, min z^t) = —1, 

X](t), x 2 (t) x 1 (f), x 2 0) 

max z 2 (t) = 4, min z 2 (t ) = 0. 

*i(f), x 2 (t ) xj 0), x 2 0) 

From the maximum and minimum values, z/t) and z 2 (t) can be represented 
by 

Zi(t) = x 1 (t)xj(t) =M l (z 1 (t))- 1 + M 1 (z 1 (t)) • ( — 1), 

^ 2 (0 = (3 + * 2 (0)*i(0 = N 1 (z 2 (t)) • 4 + N 2 (z 2 (0) • 0, 

where 

M 1 (z 1 (0)+M 2 (z 1 (0) = 1, 

N x (z 2 (t)) + N 2 (z 2 (t)) = 1. 

Therefore the membership functions can be calculated as 

Zi(f) + 1 1 - z,(f) 

M i(^i(0) = 2 ’ M 2( z i(0) = 2 ’ 

, .V * 2 (0 4-z 2 (0 

^l( Z 2(0) = ~ » ^2(^(0) = 4 • 

We name the membership functions “Positive,” “Negative,” “Big,” and 
“Small,” respectively. Then, the nonlinear system (2.10) is represented by the 
following fuzzy model. 

Model Rule 1: 

IF z,(t) is “Positive” and z 2 (t ) is “Big,” 

THEN jf (t) = /l,x(t). 

Model Rule 2: 

IF Z|(t) is “Positive” and z 2 0) is “Small,” 

THEN x(t ) =zt 2 jc(t). 

Model Rule 3: 

IF z,(t) is “Negative” and z 2 (f) is “Big,” 
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M 2 (z,(f)) M x (z x (t)) 




Fig. 2.4 Membership functions M/zjff)) and M 2 (z l (t)). 



NAzAt)) N.( zAt)) 




Fig. 2.5 Membership functions N l (z 2 (t)) and NAz 2 (t)). 



THEN x(t) = A 3 x(t). 

Model Rule 4: 

IF z x (0 is “Negative” and z 2 (t) is “Small,” 
THEN x(t) =A 4 x(t). 

Here, 




Figures 2.4 and 2.5 show the membership functions. 
The defuzzification is carried out as 

x(t) = Yj hi(z(t))AjX(t), 
i = 1 
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where 



/?i(z(0) =M 1 (z 1 (t)) XJVj(z 2 (0), 
h 2 (z(t)) =M 1 (z 1 (t )) XiV 2 (z 2 (0), 

M*( 0 ) =M 2 (zi( 0 ) X ^1(22(0), 
h 4 (z(t)) = M 2 (z 1 (z)) X N 2 (z 2 (t)). 

This fuzzy model exactly represents the nonlinear system in the region 
[- 1, 1] X [ - 1, 1] on the x r x 2 space. 

Example 3 The equations of motion for the inverted pendulum [21] are 
*1 (0 = * 2 ( 0 . 

gsin(x 1 (?)) - ambcj(t) sin(2x 1 (?))/2 - a cos(x 1 (?))u(z) 

^ 4//3 — ami cos 2 (x^t)) 

( 2 . 11 ) 



where Tj(/) denotes the angle (in radians) of the pendulum from the vertical 
and x 2 (t) is the angular velocity; g = 9.8 m /s 1 is the gravity constant, 
m is the mass of the pendulum, M is the mass of the cart, 2/ is the length 
of the pendulum, and u is the force applied to the cart (in newtons); 
a = 1 /(m + M). 

Equation (2.11) is rewritten as 



* 2 (0 



4//3 — ami cos 2 (x^t)) 

amlx 2 (t ) sin(2x 1 (t)) 



X f g sin(x 1 (t)) - 



x 2 (t) - a cos(x 1 (t))u(t) j . 

( 2 . 12 ) 



Define 



1 

Z (t) = 

1 4//3 — um/cos 2 (x 1 (t)) ’ 

z 2 (0 = sin(x!(0), 

z 3 (t) = x 2 (t) sin(2x 1 (t)), 
z 4 (t) = cos(x 1 (t)), 



where x 4 (t) e (—tt/2, tt/2 ) and x 2 (t) e [—a, a]. Note that the system is 
uncontrollable when x/t) = + 7r/2. To maintain controllability of the fuzzy 
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model, we assume that x/t) e [ — 88°, 88°]. Equation (2.12) is rewritten as 
I ami 

X 2 (t) = z 1 (t)lgz 2 (t) —z 3 (t)x 2 (t) — az 4 (t)u(t) 

As shown in Example 2, we replace z t (t) — z 4 (t) with T-S fuzzy model 
representation. Since 

1 

maxz.(t) = r = q ,, B = cos(88°), 

Xl (t) K ’ 41/3 - ami (3 2 n H V ’ 

1 



minz/t) 

x 1 (/) 



4//3 — ami 



— Qii 



z/t) can be rewritten as 



z i(0 = E £ /( z i(0)<?o 



i = 1 



(2.13) 



where 



£i( z i(0) 



z i(0 ~ <h 

<?i - <?2 



^ 2 ( z i(0) 



gt ~ z i(0 
<?1 - <?2 



The membership functions, /z ,( z ,( r)) and F/z/t)), are obtained from the 
property of E/z/t)) + E 2 (z 1 (t)) = 1. 

Figure 2.6 shows z 2 (0 = sin(xj(t)) and its local sector, where x/t) e 
( — 7 t/2, 7t/ 2). From Figure 2.6, we can find the sector [b 2 , 6, ] that consists of 
two lines b l x 1 and b 2 x v where the slopes are b 1 = 1 and b 2 = 2/7 r. 




Fig. 2.6 si n( A' // ) ) and its sector. 
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Therefore, we represent sinfx/O) as follows: 



z 2 (0 = sin(x 1 (0) = L M i(z 2 (t)) b i *i(0- 



i = 1 



From the property of membership functions + M 2 (z 2 (t )) 

can obtain the membership functions 

f z 2 (t) ~ (2 /tt) Sin~ 1 (z 2 (t)) 

(1-2/ tt) Sin - x (z 2 (0) ’ ZAt ’ ’ 
\ 1, otherwise, 

( Sin _1 (z 2 (t)) — z 2 (t) 

M 2 (z 2 {t )) = (1 - 2/ 77 )Sin _1 (z 2 (t)) ’ ^ # ° 

\ 0, otherwise. 

Next, consider z 3 (t) = x 2 (t) sin(2x 1 (t)). Since 

max z 3 (t) = a = c x and min z 3 (t) = — a = c 2 , 

X 2 (t) X](t), x 2 (t) 

we can derive in the same way as the z 3 (t) case: 



3 (0 =* 2 (0 sin(2x 1 (t)) = £A/;-(z(0)Ci, 



/= 1 



where 



z,(f) - c, Ci — z,(f) 

tVi(z 3 (t)) = * 2 (z 3 (0) = ~ 



We take the same procedure for z 4 (t ) as well. Since 

max z 4 (t) = 1 = d 1 and min z 4 (t) = (3 = d 2 , 

*i(t) *i(0 



we obtain 



Z 4 (0 = costx/O) = L^i(z(0) rf i, 

i = l 



Z 4 (0 - rf 2 _ , , “* 4 (0 

5i(z 4 (0) = , _ , » -S 2 (z 4 (0) = 



^1 ^2 



d 1 ^2 



(2.14) 
1], we 



(2.15) 



(2.16) 



where 
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From (2.13)— (2.16), we construct the following Takagi-Sugeno fuzzy model 
for the inverted pendulum: 

2 2 2 2 

= E E E E £ i( Z l(0) M ;( Z 2(0) Ar fc( Z 3 (0) 5 /( Z 4(0) 

7 = 1 y = l A: = 1 /= 1 







0 1 


Xj (0 




0 




X 




ami 


+ 




u{t ) 






8 ■ <li h j ~ 2 Qi C k 


x 2 (t) 




— a ■ qjdj 



2 2 2 2 

= E E E 'LE i (z 1 (t))Mj(z 2 (t))N k (z 3 (t))S,(z 4 (t)) 

/ = 1 y = l A = 1 /= 1 



x M.7H*(0 + ( 2 - 17 ) 

The summations in (2.17) can be aggregated as one summation: 

16 

*(0 = E hp(z(t)){A*x(t) + B*u(t)}, (2.18) 

P= 1 

where 

p = l + 2(k — 1) + 4(y — 1) + 8(/ — 1), 
fy,(z(0) =£ i( Z l(0) M ;( z 2(0) ^Ifc(^3(0) 5 /( z 4(0)» 

^J=^yw» B p= B ijkl- 

Equation (2.18) means that the fuzzy model has the following 16 rules: 
Model Rule 1: 

IF z ^t) is “Positive” and z 2 (t) is “Zero” 

and z 3 (t) is “Positive” and z 4 (t) is “Big,” 

THEN xit) = Afrit) + Bfuit). 

Model Rule 2: 

IF z jit) is “Positive” and z 2 it) is “Zero” 

and z 3 it) is “Positive” and z 4 it) is “Small,” 

THEN xit) = A%xit) + Bfrit). 

Model Rule 3: 

IF z jit) is “Positive” and z 2 it) is “Zero” 

and z 3 it) is “Negative” and z 4 it) is “Big,” 



ii(t) 

* 2 (0 



THEN xit) = Afrit) + B*uit). 
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Model Rule 4: 

IF z t it) is “Positive” and z 2 0) is “Zero” 

and z 3 it) is “Negative” and z 4 (t) is “Small,” 

THEN xit) = A*x(t) + B*u(t). 

Model Rule 5: 

IF z 3 (t) is “Positive” and z 2 (t) is “Not Zero” 

and z 3 it) is “Positive” and z 4 (t) is “Big,” 

THEN xit) = A* s xit) + 

Model Rule 6: 

IF z : (t) is “Positive” and z 2 (t ) is “Not Zero” 

and z 3 (t ) is “Positive” and z 4 (t) is “Small,” 

THEN x{t ) = A%x{t) + B*u(t). 

Model Rule 1: 

IF z ^t) is “Positive” and z 2 (0 is “Not Zero” 

and z 3 (t) is “Negative” and z 4 (t) is “Big,” 

THEN x(t) = A*x(t) + B*u(t). 

Model Rule 8: 

IF z 3 (t) is “Positive” and z 2 it) is “Not Zero” 

and z 3 it) is “Negative” and z 4 (t) is “Small,” 

THEN xit) = A%xit) + B%uit). 

Model Rule 9: 

IF z 3 it) is “Negative” and z 2 it) is “Zero” 

and z 3 it) is “Positive” and z 4 it) is “Big,” 



THEN xit) = A* g xit) + B*uit). 




CONSTRUCTION OF FUZZY MODEL 1 9 

Model Rule 10: 

IF z 3 (t) is “Negative” and z 2 (t ) is “Zero” 

and z 3 (t) is “Positive” and z 4 (t) is “Small,” 

THEN x(t) = A* 0 x(t) + B* 0 u(t). 

Model Rule 11: 

IF z x (t) is “Negative” and z 2 (t) is “Zero” 

and z 3 (t) is “Negative” and z 4 (t ) is “Big,” 

THEN x(t) = A^xit) + B£u(t). 

Model Rule 12: 

IF z : (t) is “Negative” and z 2 (t) is “Zero” 

and z 3 (t) is “Negative” and z 4 (t) is “Small,” 

THEN x(t) = A^xit) + B^uit). 

Model Rule 13: 

IF Zj(t) is “Negative” and z 2 (t ) is “Not Zero” 

and z 3 (t) is “Positive” and z 4 (t) is “Big,” 

THEN x(t) = A£x(t) + B£u(t). 

Model Rule 14: 

IF z x (t) is “Negative” and z 2 (t) is “Not Zero” 

and z 3 (t) is “Positive” and z 4 (t) is “Small,” 

THEN x(t) = A* A x(t) + B* 4 u(t). 

Model Rule 15: 

IF z : (t) is “Negative” and z 2 (t ) is “Not Zero” 

and z 3 (t) is “Negative” and z 4 (t) is “Big,” 



THEN x(t) = + B*u(t). 




to* 
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Model Rule 16: 

IF z 3 (t) is “Negative” and z 2 (t ) is “Not Zero” 

and z 3 (t) is “Negative” and z 4 (r) is “Small,” 

THEN x(t) = A*x(t) + B* 6 u{t). 

Here, z 3 (t), z 2 (t), z 3 (t) and z 4 (t) are premise variables and 



A* = A = 

A i -'Mill 



— ^1112 ~ 



A* = A = 

A 3 a U21 



A* = A = 

A 4 a 1122 



A* = A 
A 5 a 1211 



A* = A 

A 6 a 1212 



A* = A = 

A 7 a 1221 



A* = A = 

A S A 1222 



A* = A 

A 9 a 2111 



0 
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o 

g-q A 

o 

g-q\bi 

o 
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o 

g-q\ b 2 

o 

g-q\b 2 

o 
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ami 

1 
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1 
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1 

ami 

^■ q i Cl 

1 

ami 



Bf = Sim 



fi? = film 



fi* = film 



fi? = film 



fi? = film 



fi? = film 



fi? = fimi 



fi? fil222 



fi? = fimi 



0 

^ *^1^1 
0 

0 

2^2/ 2 
0 

0 

£2 *^1^1 
0 

0 
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0 

22 ' CJ -^(l ^ 
0 

'22 * 2/2^1 
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Figures 2.7-2.10 show the membership functions, that is, 



Ei(z x {t)) 



Zl(0 ~ 
<h ~ <h 



E 2 (zi(t )) 



gi ~^i(0 

< 7 l - <?2 



^i(z 2 (0) 



sin(xi(t)) - (2 /tt)z 2 (Q 
(1 - 2/7r)z 2 (t) 



M 2 (z 2 (t)) 



*i (0 - z 2 (t) 

(1 - 2/7r)z 2 (t) ’ 



^i( z 3 (0) 



* 3 (0 ~ c 2 

C 1 — C 2 



^ 2 ( Z 3 (0) 



^1 - * 3 (0 
C 1 — c 2 



* 4 (0 - d 2 

d ^ ^2 



S 2 (z 4 (t)) 



<7, - z 4 (t) 
d 2 



5i(z 4 (t)) 
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41/3 - ami 4//3 — aml/3 2 

Fig. 2.7 Membership functions E/z/t)) and E 2 (z/t)). 




Fig. 2.8 Membership functions M/z 2 (t)) and M 2 (z 2 (t)). 




Fig. 2.9 Membership functions M 1 (z 3 (f)) and N 2 (z 3 (tJ). 




Fig. 2.10 Membership functions S^z/f)) and S 2 (z 4 (t)). 



1 
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Remark 2 Prior to applying the sector nonlinearity approach, it is often a 
good practice to simplify the original nonlinear model as much as possible. 
This step is important for practical applications because it always leads to the 
reduction of the number of model rules, which reduces the effort for analysis 
and design of control systems. This aspect will be illustrated in design 
examples throughout this book. For instance, in the vehicle control described 
in Chapter 8, a two-rule fuzzy model is obtained. If we attempt to derive a 
fuzzy model without simplifying the original nonlinear model, 2 6 rules would 
be needed to exactly represent the nonlinear model. We will see in Chapter 8 
that the fuzzy controller design based on the two-rule fuzzy model performs 
well even for the original nonlinear system. 

2.2.2 Local Approximation in Fuzzy Partition Spaces 

Another approach to obtain T-S fuzzy models is the so-called local approxi- 
mation in fuzzy partition spaces. The spirit of the approach is to approximate 
nonlinear terms by judiciously chosen linear terms. This procedure leads to 
reduction of the number of model rules. For instance, the fuzzy model for 
the inverted pendulum in Example 3 has 16 rules. In comparison, in Example 
4 a 2-rule fuzzy model will be constructed using the local approximation idea. 
The number of model rules is directly related to complexity of analysis and 
design LMI conditions. This is because the number of rules for the overall 
control system is basically the combination of the model rules and control 
rules. 

Remark 3 As pointed out above, the local approximation technique leads to 
the reduction of the number of rules for fuzzy models. However, designing 
control laws based on the approximated fuzzy model may not guarantee the 
stability of the original nonlinear systems under such control laws. One of the 
approaches to alleviate the problem is to introduce robust controller design, 
described in Chapter 5. 

Example 4 Recall the inverted pendulum in Example 3. In that example, 
the constructed fuzzy model has 16 rules. In the following we attempt to 
construct a two-rule fuzzy model by local approximation in fuzzy partition 
spaces. Of course, the derived model is only an approximation to the original 
system. However, it will be shown later in this chapter that a fuzzy controller 
design based on the two-rule fuzzy model performs well when applied to the 
original nonlinear pendulum system. 

When .t | ( / ) is near zero, the nonlinear equations can be simplified as 

Xi(t)=x 2 (t), (2.19) 

g*i(Q ~ a«(Q 

2 41/3 — ami 



( 2 . 20 ) 
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When x/t) is near + tt/2, the nonlinear equations can be simplified as 



*i(0 =* 2 (0» 

2gr 1 (t) /7r - af3u{t ) 



w(0 



4//3 — ami (3 ‘ 



( 2 . 21 ) 

( 2 . 22 ) 



where j3 = cos(88°). 

Note that (2.19)-(2.22) are now linear systems. We arrive at the following 
fuzzy model based on the linear subsystems: 

Model Rule 1 



IF x ^t) is about 0, 

THEN x(t) = Afx(t) + Bfiiit). 



Model Rule 2: 

IF x f(t) is about +7r/2(|x 1 | < 7r/2), 
THEN xit) = A 2 xit) + B 2 uit). 

Here, 



A 



^2 



0 1 

— * 0> 

4//3 — ami 

0 

2g 

7 t ( 4//3 — ami 1 3 2 ) 



1 

0 ’ 






b 2 



a 

41/3 — ami 
0 

a (3 

41/3 — ami (3 2 



and /3 = cos(88°). Membership functions for Rules 1 and 2 can be simply 
defined as shown in Figure 2.11. 

Remark 4 In Example 4, the membership functions are simply defined 
using triangular types. Note that the fuzzy model is an approximated model. 
Therefore we may simply define triangular-type membership functions. On 
the other hand, in the fuzzy model in Example 3, the membership functions 
are obtained so as to exactly represent the nonlinear dynamics. 

The following remark addresses the important issue of approximating 
nonlinear systems via T-S models. 
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Rule 2 




Fig. 2.11 Membership functions of two-rule model. 



Remark 5 Section 2.2 presents the approaches to obtain a fuzzy model for a 
nonlinear system. An important and natural question arises in the construc- 
tion using local approximation in fuzzy partition spaces or simplification 
before using sector nonlinearity. One may ask, “Is it possible to approximate 
any smooth nonlinear systems with Takagi-Sugeno fuzzy models (2.1) having 
no consequent constant terms?” The answer is fortunately Yes if we consider 
the problem in C° or C 1 context. That is, the original vector field plus its 
first-order derivative can be accurately approximated. Details will be pre- 
sented in Chapter 14. 



2.3 PARALLEL DISTRIBUTED COMPENSATION 

The history of the so-called parallel distributed compensation (PDC) began 
with a model-based design procedure proposed by Kang and Sugeno (e.g., 
[16]). However, the stability of the control systems was not addressed in the 
design procedure. The design procedure was improved and the stability of 
the control systems was analyzed in [2]. The design procedure is named 
"parallel distributed compensation” in [14]. 

The PDC [2, 14, 15] offers a procedure to design a fuzzy controller from a 
given T-S fuzzy model. To realize the PDC, a controlled object (nonlinear 
system) is first represented by a T-S fuzzy model. We emphasize that many 
real systems, for example, mechanical systems and chaotic systems, can be 
and have been represented by T-S fuzzy models. 

In the PDC design, each control rule is designed from the corresponding 
rule of a T-S fuzzy model. The designed fuzzy controller shares the same 
fuzzy sets with the fuzzy model in the premise parts. For the fuzzy models 
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(2.1) and (2.2), we construct the following fuzzy controller via the PDC: 
Control Rule i: 

IF z ^t) is M a and ••• and z p (t) is M jp , 

THEN uit) = —Fjxit), i = 1,2, . . . , r. 

The fuzzy control rules have a linear controller (state feedback laws in this 
case) in the consequent parts. We can use other controllers, for example, 
output feedback controllers and dynamic output feedback controllers, instead 
of the state feedback controllers. For details, consult Chapters 12 and 13, 
which are devoted to the problem of dynamic output feedback. 

The overall fuzzy controller is represented by 

u(t) = - ' -1 r = - 'Zh i {z{t))F i x{t). (2.23) 

X>,(z(0) ,=1 

i = 1 

The fuzzy controller design is to determine the local feedback gains F l in 
the consequent parts. With PDC we have a simple and natural procedure to 
handle nonlinear control systems. Other nonlinear control techniques require 
special and rather involved knowledge. 

Remark 6 Although the fuzzy controller (2.23) is constructed using the local 
design structure, the feedback gains F t should be determined using global 
design conditions. The global design conditions are needed to guarantee the 
global stability and control performance. An interesting example will be 
presented in the next section. 

Example 5 If the controlled object is represented as the model rules shown 
in Example 1, the following control rules can be constructed via the PDC: 

Control Rule i: 

IF x(t) is M n and ••• and x(t — n + 1) is M in , 

THEN uit ) = —Fjxit), i = 1,2,..., r. 

2.4 A MOTIVATING EXAMPLE 

In this chapter, for brevity only results for discrete-time systems are pre- 
sented. The results, however, also hold for continuous-time systems subject to 
some minor modifications. 
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The open-loop system of (2.5) is 

r 

*(t+l) = Yj /i,(z(0Mi*(0- (2.24) 

i= 1 

A sufficient stability condition, derived by Tanaka and Sugeno [1, 2], for 
ensuring stability of (2.24) follows. 

THEOREM 1 [1, 2] The equilibrium of a fuzzy system (2.24) is globally asymp- 
totically stable if there exists a common positive definite matrix P such that 

A T i PA i — P < 0, i = 1, 2, . . . , r, (2.25) 

that is, a common P has to exist for all subsystems. 

This theorem reduces to the Lyapunov stability theorem for (discrete-time) 
linear systems when r = 1. 

The stability condition of Theorem 1 is derived using a quadratic function 
F(j c(t)) = x(t) T Px(t). If there exists a P > 0 such that K(j c(t)) = x(t) T Px(t) 
proves the stability of system (2.24), system (2.24) is also said to be quadrati- 
cally stable and V(x(l)) is called a quadratic Lyapunov function. Theorem 1 
thus presents a sufficient condition for the quadratic stability of system 
(2.24). 

To check the stability of fuzzy system (2.24), the lack of systematic 
procedures to find a common positive definite matrix P has long been 
recognized. Most of the time a trial-and-error type of procedure has been 
used [2, 23]. In [13] a procedure to construct a common P is given for 
second-order fuzzy systems, that is, the dimension of state n = 2. We first 
pointed out in [14, 15, 24] that the common P problem can be solved 
efficiently via convex optimization techniques for LMIs [18]. To do this, a 
very important observation is that the stability condition of Theorem 1 is 
expressed in LMIs. To check stability, we need to find a common P or 
determine that no such P exists. This is an LMI problem. See Section 2.5.2 
for details on LMIs and the related LMI approach to stability analysis and 
design of fuzzy control systems. Numerically the LMI problems can be solved 
very efficiently by means of some of the most powerful tools available to date 
in the mathematical programming literature. For instance, the recently 
developed interior-point methods [19] are extremely efficient in practice. 

A question naturally arises of whether system (2.24) is stable if all its 
subsystems are stable, that is, all A- s are stable. The answer is no in general, 
as illustrated by the following example. 

Example 6 Consider the following fuzzy system: 

Rule 1: 



IF x 2 (t ) is My (e.g., Small), 
THEN x{t + 1) =A lX (t). 
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Fig. 2.12 Membership functions of Example 6. 



Rule 2: 

IF x 2 (t) is M 2 (e.g., Big), 

THEN x(t + 1) =A 2 xU). 
Here, x(t) = [x^t) x 2 (t)] T and 



1 


-0.5' 


, An = 


-l 


-0.5' 


1 


0 




l 


0 



Figure 2.12 shows the membership functions of M x and M 2 . Since A , and 
A 2 are stable, the linear subsystems are stable. However, for some initial 
conditions the fuzzy system can be unstable, as shown in Figure 2.13 for the 
initial condition x = [0.90 — 0.70] r . It should be noted that the linearization 
of the fuzzy system around 0 is stable (which implies that the fuzzy system is 
locally stable). Obviously there does not exist a common P > 0 since the 
fuzzy system is unstable. This can be shown analytically. Moreover this can 
also be shown numerically by convex optimization algorithms involving LMIs. 

Still an interesting question is for what initial conditions the fuzzy system 
is stable (or unstable). This is determined by studying the basin of attraction 
of the origin. 1 

Figure 2.14(a) shows the basin of attraction for the case of a = 1. The 
black area indicates regions of instability (horizontal axis is x , ). It is also of 
interest to consider how the basin of attraction changes as the membership 
functions vary, for instance, how the basin of attraction would change as a 
varies for this example. Figures (b), (c), and (d) show the basin of attraction 

1 Sugeno mentioned this point in his plenary talk titled “Fuzzy Control: Principles, Practice, and 
Perspectives” at 1992 IEEE International Conference on Fuzzy Systems, March 9, 1992. 
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Fig. 2.13 Response of Example 6 (a = 1). 



/ 

/ 




- 

| 



for various values of a. It can be seen that as a decreases (increases) from 1, 
the basin of attraction becomes smaller (larger). Therefore, the basin of 
attraction for the fuzzy system could be membership function dependent. In 
the example, when a = °°, the fuzzy system becomes 

A-i + A . 7 

x (t + 1 ) = 2 — ~ x ^’ 

which is linear and globally asymptotically stable. 

For this example, an interesting interpretation can be given for the 
dependence of basin of attraction on membership functions. As a increases 
(decreases), the inference process tends to be “fuzzier” (“crisper”). Hence a 
fuzzier decision leads to a larger basin of attraction while a crisper decision 
leads to a smaller basin of attraction. 

As illustrated by the example, we have to take stability into consideration 
when selecting rules and membership functions. How to systematically select 
rules and membership functions to satisfy prescribed stability properties is an 
interesting topic. In the next section, we consider the control design problems 
via parallel distributed compensations. 



2.5 ORIGIN OF THE LMI-BASED DESIGN APPROACH 

This section gives the origin of the control design approach, which forms the 
core subject of this book, that is, the LMI-based design approach. The 
objective here is to illustrate the basic ideas [24] of stability analysis and 
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-2.5 -x, - 2.5 -2.5 -x 2 - 2.5 

(a) a = 1 




-2.5 -x, - 2.5 -2.5 -x 2 - 2.5 

(c) a = 0.25 




-2.5 -x, - 2.5 -2.5 — x 2 - 2.5 

( b ) a = 0.5 




I I 

-2.5 — x ! - 2.5 -2.5 -x 2 - 2.5 

(d) a = 2.0 



Fig. 2.14 Basin of attraction for Example 6. 



stable fuzzy controller design via LMIs. The details will be presented in 
Chapter 3. 

2.5.1 Stable Controller Design via Iterative Procedure 

The PDC fuzzy controller is 



u(t ) = - h i (z( t ))F i x(t). (2.26) 

i= 1 

Note that the controller (2.26) is nonlinear in general. 

Substituting (2.26) into (2.5), we obtain 



*(t+i)= E E / hO(0) /! jO(0)K' - ( 2 - 27 ) 

i= 1 7=1 

Applying Theorem 1, we have the following sufficient condition for 
(quadratic) stability. 
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THEOREM 2 The equilibrium of a fuzzy control system (2.27) is globally 
asymptotically stable if there exists a common positive definite matrix P such that 



{A, - B i F j } T P{A i - BjFj) -P< 0 



(2.28) 



for hj(z(t )) • hj(z(t )) ¥= 0, \/t, i, j = 1, 2, . . . , r. 
Note that system (2.27) can also be written as 



x{t + 1) 



X>,(z(0)/?,(z(0)M, - B t F t }x(t) 



i= 1 



+ 2 £ 'Lhi{z{t))h j {z{t))G ij x{t) 

i = 1 i<j 



(2.29) 



where 



[A, - B,F + {/l, - «,/•,} 

Gy = — - 1 J , i <j S.t. A,, n h j * </>• 



Therefore we have the following sufficient condition. 



THEOREM 3 77ie equilibrium of a fizzy control system (2.27) is globally 
asymptotically stable if there exists a common positive definite matrix P such that 
the following two conditions are satisfied : 

{A, - B i F l ) T P{A l - BiF t ) - P < 0, i= 1,2,..., r (2.30) 
GfjPGjj - P < 0, i <j < r s.t. h t n h j ¥= f. (2.31) 

For the meaning of the notation i < j < r s.t. h t n hj ¥= <h, see 
Chapter 1. 



Remark 7 The conditions of Theorem 3 are more relaxed than those of 
Theorem 2. 



The control design problem is to select F t (i = 1, 2, ... , r) such that 
conditions (2.30) and (2.31) in Theorem 3 are satisfied. Using the notation of 
quadratic stability, we can also formulate the control design problem as to 
find F’s such that the closed-loop system (2.27) is quadratically stable. 
If there exist such F’s, the system (2.5) is also said to be quadratically 
stabilizable via PDC design. 

In this chapter, we first design a controller for each rule and check 
whether the stability conditions are satisfied. Recall we can use LMI convex 
programming techniques to solve this stability analysis problem. If the stabil- 
ity conditions are not satisfied, we have to repeat the procedure. Consult 
Section 2.5.2 on how LMIs can be used to directly solve the control design 
problem. 
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Next consider the common B matrix case, that is, B t = B(i = 1, 2 , . . . , r). 
In this case, Theorem 3 reduces to: 

THEOREM 4 When B t = B, i = 1, . . . , r, the equilibrium of the fuzzy control 
system (2.27) is globally asymptotically stable if there exists a common positive 
definite matrix P such that 

[A i - B i F l ) T P{A i - B,F t ) - P < 0, f = 1,2,..., r. (2.32) 
Furthermore, for the common B case, if we can choose F t such that 

A t - BFj = G, (2.33) 

where G is a Hurwitz matrix, then the system (2.27) becomes a linear system 

x(t + 1) = Gx(t). 

This is a global linearization result. We remark that a common G might not 
always be possible even if (A t , Bf) are controllable. 

Remark 8 As shown in Theorem 4, the stability conditions are simplified in 
the common B matrix case. The same feature will be observed in all the 
chapters. 

Let us look at some examples. 

Example 7 Consider the following fuzzy system: 

Model Rule 1: 

IF x 2 (t) is Mj, 

THEN x(t + 1) = A lX (t) + Bu(t). 

Model Rule 2: 

IF x 2 (t) is M 2 , 

THEN x{t + 1) = A 2 x(t) + Bu(t). 

Here, A 1 , A 2 are the same as in Example 6 and 
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Employ the PDC controller (2.26) and choose the closed-loop eigenvalues 
to be [0.5 0.35]. We obtain 

F 1 = [0.15 -0.3250], 

F 2 = [- 1.85 -0.3250], 



and 



A | — BF — A 2 — BF 2 — G — 



0.85 

1 



-0.1750 

0 



The closed loop becomes 

x(t + 1) = Gx(t ), 

which is stable since G is stable. 

Next we consider the more general case. 

Example 8 

Consider the following fuzzy system: 

Model Rule 1: 

IF x 2 (t) is M 1; 

THEN x(t + 1) = Aj^xit) + Bpiit). 
Model Rule 2: 

IF x 2 (t) is M 2 , 

THEN x(t + 1) = Atx( 0 + B 2 u(t). 
Here, A x , A 2 are the same as in Example 6 and 



1 


, B , = 


' - 2 ' 


1 


? 2 


1 



The membership functions of Example 6 (a = 1) are used in the simulation. 
Again choose the closed-loop eigenvalues to be [0.5 , 0.35]. We have 

F l = [0.65 -0.5], 

F 2 = [0.87 -0.11], 
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and 



0.35 0 ' 


, A- 


7 — ^7 F 1 = 


'0.74 


-0.72' 


0.35 0.5. 




2 2 2 


0.13 


0.11 


0.2150 - 


0.9450' 








0.2400 


0.3050 









Note that G u is stable. 

The PDC controller is given as follows: 

Control Rule 1: 



IF x 2 (t) is M 1; 

THEN u(t) = -F lX (t). 

Control Rule 2: 

IF x 2 (t) is M 2 , 

THEN u(t) = —F 2 x(t). 

It can be easily shown that if we choose the positive definite matrix P to 
be 

[ 1.1810 -0.0614' 

[-0.0614 2.3044 ]’ 

the stability conditions (2.30) and (2.31) are satisfied. In other words, the 
closed-loop fuzzy control system which consists of the fuzzy model and the 
PDC controller is globally asymptotically stable. The P is obtained by 
utilizing an LMI optimization algorithm. Figure 2.15 illustrates the behavior 
of the fuzzy control system for the same initial condition of Figure 2.13. 

In the next section, we present an introduction to LMIs as well as the LMI 
approach to stability analysis and design of fuzzy control systems. 



2.5.2 Stable Controller Design via Linear Matrix Inequalities 

Recently a class of numerical optimization problems called linear matrix 
inequality (LMI) problems has received significant attention [18]. These 
optimization problems can be solved in polynomial time and hence are 
tractable, at least in a theoretical sense. The recently developed interior-point 
methods [19] for these problems have been found to be extremely efficient in 
practice. For systems and control, the importance of LMI optimization stems 
from the fact that a wide variety of system and control problems can be 
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Fig. 2.15 Response of Example 8. 
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recast as LMI problems [18]. Except for a few special cases these problems 
do not have analytical solutions. However, the main point is that through the 
LMI framework they can be efficiently solved numerically in all cases. 
Therefore recasting a control problem as an LMI problem is equivalent to 
finding a “solution” to the original problem. 

DEFINITION 1 [18] An LMI is a matrix inequality of the form 

m 

F(x)=F 0 + 'Zx i F i > 0, (2.34) 

i = 1 

where x T = (x t , x 2 , . . . , x m ) is the variable and the symmetric matrices F t = 
Ff e R” x ", i = 0, . . . , m, are given. The inequality symbol > 0 means that 
Fix) is positive definite. 

The LMI (2.34) is a convex constraint on x, that is, the set {x|L(jr) > 0} is 
convex. The LMI (2.34) can represent a wide variety of convex constraints on 
x. In particular, linear inequalities, convex quadratic inequalities, matrix 
norm inequalities, and constraints that arise in control theory, such as 
Lyapunov and convex quadratic matrix inequalities, can all be cast in the 
form of an LMI. Multiple LMIs > 0, i = 1 ,..., p, can be expressed as a 
single LMI diag(L (1) , . . . , F (p) ) > 0. 

Very often in the LMIs the variables are matrices, for example, the 
Lyapunov inequality 



A t PA - P < 0, 



(2.35) 
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where A e [R" x " is given and P = P T is the variable. In this case the LMI 
will not be written explicitly in the form F(x ) > 0. In addition to saving 
notation, this may lead to more efficient computation [18]. Of course, the 
inequality (2.35) can be readily put in the form (2.34): take /’„ = 0, F, = 
—A r P i A + P t , where P v . . . , P m are a basis for symmetric n X n matrices. 

LMI problems [18] Given an LMI F{x) > 0, the LMI problem is to find 
x lcas such that F( x feas ) > 0 or determine that the LMI is infeasible. This is a 
convex feasibility problem. 

As an example, the simultaneous Lyapunov stability condition in Theorem 
1 is exactly an LMI problem: Given A { e [R” x ", i = 1, . . . , r, we need to find 
P satisfying the LMI 

P > 0, A]PA i -P< 0, i=l,2, ...,r, 

or determine that no such P exists. 

The LMI problems are tractable from both theoretical and practical 
viewpoints: They can be solved in polynomial time, and they can be solved in 
practice very efficiently by means of some of the most powerful tools 
available to date in the mathematical programming literature (e.g., the 
recently developed interior-point methods [19]). 

The stability conditions encountered in this book are expressed in the 
form of LMIs. This recasting is significant in the sense that efficient convex 
optimization algorithms can be used for stability analysis and control design 
problems. The recasting therefore constitutes solutions to the stability analy- 
sis and control design problems in the framework of the Takagi-Sugeno fuzzy 
model and PDC design. 

The design procedure presented in the previous section involves an 
iterative process. For each rule a controller is designed based on considera- 
tion of local performance only. Then an LMI-based stability analysis is 
carried out to check whether the stability conditions are satisfied. In the case 
that the stability conditions are not satisfied, the controller for each rule will 
be redesigned. The iterative design procedure has been very effective in our 
experience. However, from the standpoint of control design, it is more 
desirable to be able to directly design a control that ensures the stability of 
the closed-loop system. This is referred as the control problem in the 
framework of the Takagi-Sugeno fuzzy model and PDC design. We claim 
that the control problem can be recast (hence solved) using the LMI 
approach. Here we only briefly state the ideas of the LMI approach to the 
control design problem. We show a simple case (r = 1), that is, the linear 
case, below. Fuzzy control case will be presented in Chapter 3. 

Consider the case r = 1, that is, there is only one IF-THEN rule; (2.5) 
becomes a linear time-invariant system, 

x(t + 1) = Ax(t) + Bu(t). 



(2.36) 
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For a given control gain F, using standard stability theory for linear time-in- 
variant systems or Theorem 2, the system (2.36) is (quadratically) stable if 
there exists P > 0 such that 

[A - BF} t P{A - BF] - P < 0 . (2.37) 

The control design problem is to find a state feedback gain F such that 
the closed-loop system is (quadratically) stable. If such a gain F exists, the 
system is said to be quadratically stabilizable (via linear state feedback). This 
quadratic stabilizability problem can be recast as an LMI problem. 

The condition (2.37) is not jointly convex in F and P. Now multiplying the 
inequality on the left and right by P 1 , and defining a new variable X = P 1 , 
we may rewrite (2.37) as 

X{A - BF} T X~ l {A - BF}X-X< 0 . (2.38) 

Define M = FX so that for X > 0 we have F = MX~ l . Substituting into 
(2.38) yields 



X - {AX - BM} T X l {AX - BM) > 0 . (2.39) 

This nonlinear (convex) inequality can now be converted to LMI form using 
Schur complements [18]. The resulting LMI is 



X 

{AX- BM) 



{AX- BM) t 
X 



(2.40) 



in X and M. Thus the system (2.36) is quadratically stabilizable if there exist 
X > 0 and M such that the LMI (2.40) holds. The state feedback gain is 
F = MX 1 . 

We can easily extend the LMI-based control design approach to multiple- 
rule (r > 1) cases of the Takagi-Sugeno fuzzy models. For instance, the 
quadratic stabilizability of the Takagi-Sugeno fuzzy models via a linear state 
feedback can be cast as the following LMI problem in X and M: 



X > 0 , 



X {AjX - BjM) 

{A'X-B'M) X 



>0, i = 1, 2, . . . , r. 



with the state feedback gain F = MX 

The LMI-based control design approach has also been developed for the 
control of Takagi-Sugeno fuzzy models via PDC design. For more details, see 
Chapter 3. 

Some important remarks are in order. 
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Remark 9 The stability conditions presented in this book not only guaran- 
tee stability of fuzzy models and fuzzy control systems, they also guarantee 
stability for related uncertain linear time-varying [linear differential inclusion 
(LDI)] systems and nonlinear systems satisfying some global or local sector 
conditions. Thus a controller that works well with the fuzzy model is likely to 
work well when applied to the real system. This point is clearly demonstrated 
by the application in the next section. The theoretical details, however, will 
be discussed in other chapters. 

Remark 10 The stability analysis and control design results presented in this 
section hold for continuous-time systems as well. Instead of using the 
Lyapunov inequality for discrete-time systems, we should use the Lyapunov 
inequality for continuous-time systems, 

A t P + PA < 0 . 

In the next section, we apply the PDC approach to a continuous-time 
system. 



2.6 APPLICATION: INVERTED PENDULUM ON A CART 



To illustrate the PDC approach, consider the problem of balancing and 
swing-up of an inverted pendulum on a cart. Recall the equations of motion 
for the pendulum [21]: 



*i(0 = x 2( 0 



MO 



gsin(x 1 (t)) - amlxj(t)sin(2x 1 (t))/2 - acos(x/t))u(t) 
41/3 — ami cos 2 (x 1 (t)) 



(2.41) 



where x/t) denotes the angle (in radians) of the pendulum from the vertical 
and x 2 (t) is the angular velocity; g = 9.8 m/s 2 is the gravity constant, 
m is the mass of the pendulum, M is the mass of the cart, 2/ is the length of 
the pendulum, and u is the force applied to the cart (in newtons); and 
a = 1 /(m + M). We choose m = 2.0 kg, M = 8.0 kg, 21 = 1.0 m in the 
simulations [20]. 



2.6.1 Two-Rule Modeling and Control 

The control objective of this subsection is to balance the inverted pendulum 
for the approximate range XjE(— 7t/2, 7t/2). In order to use the PDC 
approach, we must have a fuzzy model which represents the dynamics of the 
nonlinear plant. Therefore we first represent the system (2.41) by a Takagi- 
Sugeno fuzzy model. To minimize the design effort and complexity, we try to 




APPLICATION: INVERTED PENDULUM ON A CART 39 



use as few rules as possible. Notice that when x 1 = + 7r/2, the system is 
uncontrollable. Hence, as shown in Example 4, we approximate the system by 
the following two-rule fuzzy model: 

Rule 1: 

IF x/t) is about 0, 

THEN x(t) = A 1 x(t) + Byiiit). 



Rule 2: 

IF x ^t) is about + 7r/2 (|x : | < 7 t/2), 
THEN x(t) = A 2 x(t) + B 2 u(t). 



Here, 



A 



i 



0 

8 

41/3 — ami 



1 

0 ’ 



B i = 



0 

a 

41/3 — ami 



A 2 



0 

2g 

7 t ( 4//3 — ami f3 2 ) 



1 

0 ’ 



b 2 



0 

a (3 

41/3 — ami [3 2 



and (3 = cos(88°). 

Membership functions for Rules 1 and 2 are shown in Figure 2.f6. 



Rule 2 




Fig. 2.16 Membership functions of two-rule model. 
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Choose the closed-loop eigenvalues [ — 2, —2] for A 1 — B l F l and A 2 — 
B 2 F 2 . We have 

F 1 = [-120.6667 -22.6667], 

F 2 = [-2551.6 -764.0], 



It follows that 



and 



A i B l F 1 A 2 B 2 F 2 G 



0 1 
-4 -4 



^12 — 



0 1 
-212.1325 -67.4675 



Note that G n is Hurwitz. 

Using an LMI optimization algorithm, we obtain 



3.6250 0.6250 

0.6250 0.2812 ' 



(2.42) 



It can be easily shown that the following stability conditions are satisfied: 

{A t - BjFj} 7 P + P{A, - B ,F,} <0, i = 1,2, (2.43) 

G\ 2 P + PG n < 0. (2.44) 

The resulting PDC control law is as follows: 



Rule 1: 



IF x ^t) is about 0, 

THENu(t) = -Ffx(t). 



Rule 2: 

IF x f(t) is about + 7r/2 (| JCj | < rr/2), 
THENuit) = -F 2 xit). 



That is, 

u(t) = -hf(Xf{t))FfXit) - h 2 (Xfit))F 2 xit). (2.45) 

The membership values of Rules 1 and 2 are hf and h 2 , respectively 
ihf + h 2=1). This (nonlinear) control law guarantees the stability of the 
fuzzy control system (fuzzy model + PDC control). To assess the effective- 
ness of the PDC controller, we apply the controller to the original system 
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Fig. 2.17 Angle response using linear and two-rule fuzzy control. 



(2.41). As pointed out in Remark 3, we may design a robust fuzzy controller 
(see Chapter 5 for the details) that can compensate the approximation error. 

Simulations indicate the control law can balance the pendulum for initial 
conditions x 1 e [ — 88°, 88°] (x 2 = 0). In contrast, the linear control alone 
u = —F 1 x fails to balance the pendulum for initial angles \x 1 \ > 45°. Figure 
2.17 shows the response of the pendulum system using linear and fuzzy PDC 
controls for initial conditions x 1 = 15°, 30°, 45°, and x 2 = 0. The solid lines 
indicate responses with the fuzzy controller. The dotted lines show those with 
the linear controller. Figure 2.18 illustrates the closed-loop behavior of the 
system with the fuzzy controller for initial conditions x t = 65°, 75°, 85°, and 
*2 = 0. 

We remark that given the nonlinear plant (2.41) nonlinear control laws can 
be designed to balance the pendulum for initial angles x t e ( — tt/2, tt/2). 
However, such control laws often tend to be quite involved. For example, one 
such control law is [20] u = k(x v x 2 ), where 



8 4 le x e 2 

k(x 1 ,x 2 ) = tan(x x ) — 1 ln[sec(x 1 ) + tan(xj)] 



+ e l e 2 ml sin(x : ) 



( e i + e 2 ) x 2 



41 

— sec(xj) - ami cos(xj) 



(2.46) 



and e v e 2 are the specified closed-loop eigenvalues. 





42 TAKAGI-SUGENO FUZZY MODEL AND PARALLEL DISTRIBUTED COMPENSATION 




In contrast, the PDC design is intuitive and simple. The resulting con- 
troller is simple as well. 

2.6.2 Four-Rule Modeling and Control 

Suppose the pendulum on the cart system is built in such a way that the work 
space of the pendulum is the full circle [— 77, 7 t ]. In this subsection, we 
extend the results to the range of Xj e [ — 77 7r] except for a thin strip near 
+ 7 t/2 . Balancing the pendulum for the angle range of 7 t/2 < | x 1 1 < 7 r is 
referred to as swing-up control of the pendulum. Recall that for x l = + 7 t/2 
the system is uncontrollable. We add two more rules (Rules 3 and 4) to the 
fuzzy model. 

Rule 1: 

IF x/f ) is about 0, 

THEN x(t) = A x x(t) + B^uit). 



Rule 2: 

IF Xj(t) is about +77/2 (|x 1 (t) | < 77/2), 



THEN x(t) = A 2 x(t) + B 2 u(t). 
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Rule 3: 

IF x/t) is about ±tt/2 (\x 1 (t)\> tt/2), 
THEN x(t) = A 3 x(t ) + B 3 u(t). 

Rule 4: 

IF x/t) is about tt, 



THEN x ( t ) = A 4 x(t) + B 4 u(t). 



Here A l , B 4 , A 2 , B 2 are the same as above and 



^3 



0 

2g 

7 t ( 4//3 — ami (3 2 ) 



1 

0 ’ 



B 3 



a 4 



0 1 

0 oj’ 



b 4 



0 

a (3 

4//3 — ami 1 8 2 

0 

a 

41/3 — ami 



The membership functions of this four-rule fuzzy model are shown in Figure 
2.19. 

Again choose the closed-loop eigenvalues [-2,-2] for A, — B 3 F 3 and 
A 4 - B 4 F 4 . We have 



F 3 = [2551.6 764.0], 

F 4 = [22.6667 22.6667], 



Rule 4 Rule 4 




Fig. 2.19 Membership functions of four-rule model. 




44 



TAKAGI-SUGENO FUZZY MODEL AND PARALLEL DISTRIBUTED COMPENSATION 



It follows that 

A 3 — B 3 F 3 = A 4 — B 4 F 4 = G 

and 

r = [ 0 1 

34 [-220.5230 - 67.4675 ' 

Note that G 34 is Hurwitz. 

It can be shown that the P of (2.42) satisfies the additional stability 
conditions 

{A t - B^fp + P{Aj - B,F,} < 0 , i = 3, 4, (2.47) 

G T 34 P + PG 34 < 0 . (2.48) 

There is no overlap between membership values h t and h 3 , h 4 and h 4 , h 2 
and h 3 , and h 2 and h 4 . Hence only G u and G 34 are needed in stability 
check. 

The PDC controller is given as follows: 

Rule 1: 

IF x 4 (t) is about 0, 

THENu(t) = —F l x(t). 

Rule 2: 

IF x/t) is about ±tt/2 (\x 1 (t)\< tt/2), 

THEN u(t) = —F 2 x(t). 

Rule 3: 

IF x/t) is about ±tt/2 (\x 1 (t)\> tt/2), 

THENu(t) = —F 3 x(t). 

Rule 4: 

IF x/t) is about t t, 

THEN u(t) = -F 4 x ( t ). 



That is, 



u(t) = —h 1 (x 1 (t))F 1 x(t) - h 2 (x x (t))F 2 x{t) 

- / * 3 (^i( 0) F 3- c (0 - >u(xi(t))F 4 x(t). 



(2.49) 
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This control law guarantees stability of the fuzzy control system (four-rule 
fuzzy model + PDC control). This controller is applied to the original system 
(2.41) for evaluation of its performance. Simulation results demonstrate that 
the controller (2.49) is able to balance the pendulum for all initial angles 
except when x^t) is in a thin strip 88° < |x x (t) | < 94°. The size of this thin 
strip can be reduced by adding more rules to the model and controller. 
Figure 2.20 illustrates the response of the closed-loop system for initial 
conditions x 1 = 125°, 145°, 165°, 180° and x 2 = 0. 

Note that the nonlinear controller (2.46) does not apply for 7 t/2 

< \X 1 (t) | < 77. 

Some comparisons between the linear, nonlinear, and fuzzy control de- 
signs are summarized loosely in Table 2.1. 

To test the robustness of this controller, the following simulations are 
conducted: (1) m is changed from 2.0 to 4.0 kg, (2) M is changed from 8.0 to 
4.0 kg, and (3) 21 is changed from 1.0 to 0.5 m. For each case, we simulate 



TABLE 2.1 Comparisons of Different Control Designs 





Work range 


Simple? 


Stability 


Linear 


(—77/4 77/4) 


Yes 


Local 


Nonlinear 


(-77/2 77/2) 


No 


Nonlocal 


Fuzzy PDC 


[ — 77 77] 


Yes 


Nonlocal 
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Time (sec) 



Fig. 2.23 Closed-loop angle response with / changed. 



the closed-loop system for the following initial conditions x 1 = 45°, 85°, 145°, 
180° and x 2 (t ) = 0. The results are shown in Figures 2.21, 2.22, and 2.23, 
respectively, for cases 1, 2, and 3. 

Robustness is not considered in this design. Robust fuzzy control design in 
Chapter 5 is applicable to this system. 
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LMI CONTROL PERFORMANCE 
CONDITIONS AND DESIGNS 



The preceding chapter introduced the concept and basic procedure of 
parallel distributed compensation and LMI-based designs. The goal of this 
chapter is to present the details of analysis and design via LMIs. This chapter 
forms a basic and important component of this book. To this end, it will be 
shown that various kinds of control performance specifications can be repre- 
sented in terms of LMIs. The control performance specifications include 
stability conditions, relaxed stability conditions, decay rate conditions, con- 
strains on control input and output, and disturbance rejection for both 
continuous and discrete fuzzy control systems [1-3]. Other more advanced 
control performance considerations utilizing LMI conditions will be pre- 
sented in later chapters. 



3.1 STABILITY CONDITIONS 

In the 1990’s, the issue of stability of fuzzy control systems has been 
investigated extensively in the framework of nonlinear system stability [1-18]. 
Today, there exist a large number of papers on stability analysis of fuzzy 
control in the literature. This section discusses some basic results on the 
stability of fuzzy control systems. 

In the following, Theorems 5 and 6 deal with stability conditions for the 
open-loop systems. Theorem 5 can be readily obtained via Lyapunov stability 
theory. The proof of Theorem 6 is given in [4, 7]. 
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THEOREM 5 [CFS] The equilibrium of the continuous fuzzy system (2.3) with 
u(t ) = 0 is globally asymptotically stable if there exists a common positive definite 
matrix P such that 



A]P + PAi < 0, i = 1, 2, . . . , r, (3.1) 

that is, a common P has to exist for all subsystems. 

THEOREM 6 [DFS] The equilibrium of the discrete fuzzy system (2.5) with 
u(t ) = 0 is globally asymptotically stable if there exists a common positive definite 
matrix P such that 



A]PA i -P< 0, i = 1,2, . . . , r, (3.2) 

that is, a common P has to exist for all subsystems. 

Next, let us consider the stability of the closed-loop system. By substituting 
(2.23) into (2.3) and (2.5), we obtain (3.3) and (3.4), respectively. 

CFS 

*(0 = E E fci(z(0)fy(z(0)Ui - B iFj}x(t). (3.3) 

i = 1 _/ = l 

DFS 

x(t+l) = E EA/(z(0) A ;(z(0)M«--®/^}*(0- ( 3 - 4 ) 

i = 1 7=1 



Denote 



Gij ~ - /F /';. 



Equations (3.3) and (3.4) can be rewritten as (3.5) and (3.6), respectively. 
CFS 

x(t) = I>;(z(0)M z (0) G ;;*(0 



1 = 1 



+ 2 E Lhi(z(t))hj(z(t)) 

i = 1 i<j 



G u + G n 



x(t). (3.5) 



DFS 



x(t + 1 ) = ’Lh i (z(t))hfz(t))G ii x(t) 

i = 1 



+ 2 E E*i(z(0) A ;(z(0) 

» = 1 ><j 



G a + G ji 



x(t). (3.6) 
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By applying the stability conditions for the open-loop system (Theorems 5 
and 6) to (3.5) and (3.6), we can derive stability conditions for the CFS and 
the DFS, respectively. 



THEOREM 7 [CFS] The equilibrium of the continuous fuzzy control system 
described by (3.5) is globally asymptotically stable if there exists a common 
positive definite matrix P such that 



GhP + PG U < 0 , 



(3.7) 



Gii + G ji 



Gjj + G n 

P + P | — — — - | < 0, 



i < j s.t. lij n hj # 4>- (3-8) 



Proof. It follows directly from Theorem 5. 



For the explanation of the notation i < j s.t. h t n hj # <l>. refer to 
Chapter 1. 



THEOREM 8 [DFS] The equilibrium of the discrete fizzy control system 
described by (3.6) is globally asymptotically stable if there exists a common 
positive definite matrix P such that 



GfjPGjj - P < 0 , 



(3.9) 



Gjj + Gjj \ j Gjj + Gjj 

1 1 1 P\ 1 . - 1 - P< 0, 



i < j s.t. hj n hj # fi. (3.10) 



Proof. It follows directly from Theorem 6. 

The fuzzy control design problem is to determine FI s ( j = 1,2,..., r) 
which satisfy the conditions of Theorem 7 or 8 with a common positive 
definite matrix P. 

Consider the common B matrix case, that is, B x = B 2 = ••• = B r . In this 
case, the stability conditions of Theorems 7 and 8 can be simplified as 
follows. 

COROLLARY 1 Assume that /f, = B 2 = ••• = B r . The equilibrium of the 
fuzzy control system (3.5) is globally asymptotically stable if there exists a 
common positive definite matrix P satisfying (3.7). 

COROLLARY 2 Assume that /I, = B 2 = ••• = B r . The equilibrium of the 
fuzzy control system (3.6) is globally asymptotically stable if there exists a 
common positive definite matrix P satisfying (3.9). 
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In other words, the corollaries state that in the common B case, Gj t P + 
PG U < 0 implies 



G u + G n 



p + p 



G ij + G n 



< o 



and GjjPG n - P < 0 implies 



G ij + G ji 



G U + G ji 



- P < 0 



To check stability of the fuzzy control system, it has long been considered 
difficult to find a common positive definite matrix P satisfying the conditions 
of Theorems 5-8. A trial-and-error type of procedure was first used [4, 7, 9]. 
In [19], a procedure to construct a common P is given for second-order fuzzy 
systems, that is, the dimension of the state is 2. It was first stated in 
[11, 12, 17] that the common P problem for fuzzy controller design can be 
solved numerically, that is, the stability conditions of Theorems 5-8 can be 
expressed in LMIs. For example, to check the stability conditions of Theorem 
7, we need to find P satisfying the LMIs 



P > 0, 

( G a + G n 



GjjP + PG h < 0, 



T 

P + P 



G,; + G,: 



< o, 



i < j s.t. hj n hj # 4>, 



or determine that no such P exists. This is a convex feasibility problem. As 
shown in Chapter 2, this feasibility problem can be numerically solved very 
efficiently by means of the most powerful tools available to date in the 
mathematical programming literature. 



3.2 RELAXED STABILITY CONDITIONS 

We have shown that the stability analysis of the fuzzy control system is 
reduced to a problem of finding a common P. If r, that is the number of 
IF-THEN rules, is large, it might be difficult to find a common P satisfying 
the conditions of Theorem 7 (or Theorem 8). This section presents new 
stability conditions by relaxing the conditions of Theorems 7 and 8. Theorems 
9 and 10 provide relaxed stability conditions [1-3]. First, we need the 
following corollaries to prove Theorems 9 and 10. 

COROLLARY 3 

Z>t(z( 0) - — T E E2/l i (z(0)/ly(z(0) >o, 

i = 1 r ^ i = l i<j 




